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The stability of self-similar power-law disks to non-axially symmetrical dis- 
turbances of zero frequency (u> = ) is discussed. It is established that marginally 
unstable modes either have uj = or consist of continua with modes of all fre- 
quencies becoming unstable together. A careful study of all modes close to u = 
yields a remarkable conclusion. 
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1. Introduction. 

Following early developments by Brodetsky & Snow, Kalnajs has given a 
pretty development of the theory of gravitational potential of equiangular spiral 
waves. Donner, in a fine but unpublished Cambridge thesis ( supervised by DLB), 
gave an account of viscous spiral waves in both the linear and non-linear regimes. 
The relationship of the vanishing viscosity limit of these waves to the stability 
problem of spiral structure was not completely elucidated. Zang and Toomre have 
explored the stability of the V = constant Mestel disk in stellar dynamics. Zang 
found m = 1 more unstable than the m = 2 modes. Toomre used the stable modes 
to illustrate very beautifully the effects of the swing amplifier. 

Spruit in a more recent paper that parallels Donner's analysis has discussed 
the non-linear equiangular spiral viscous accretion disk. 

We were led to take a new look at the inviscid problem because the dimen- 
sional argument of section 3 showed that there is no characteristic frequency of 
a 'power-law' disk. It was then unclear how such a disk could pick out a charac- 
teristic frequency for its marginally stable mode. This paper finds all the u> = 
modes and resolves this question. 
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2. Self-Similarity. 



A flat disk is said to be self-similar if any physical quantity q defining the 
configuration at the point R, is related to that quantity at other points by 

q(XR,<P + f(X)) = A 1 (X)q(R,<j ) ). (2.1) 

So the physical quantity at a point A times further out and rotated by /(A) is 
a multiple of that quantity at the original point. When A = 1 we ask that the 
transformed point be coincident with the original one so we take /(l) = and 
Ai(l) = 1. Furthermore, while relationship (1) shall be true of all the physical 
quantities defining the configuration and for all values of A, nevertheless the func- 
tion Ai(X) will depend on the physical dimensions mass, length, time, etc. of q 
and be such that Ai(X) = 1 if q is dimensionless . /(A) does not depend on q. 



Differentiating (2.1) with respect to A and then putting A = 1, we deduce 

, 9q dq_ 

'dR d(f> 



where a = f(l) and b = A[(l). 

Equation (2.2) may be rewritten 



R dR +a ^ = °> (2 - 3) 

where 

Q = lnq-blnR = ln(qR- b ). (2.4) 

The general solution of (2.3) is 

Q = S{(j)-alnR) , (2.5) 

where S is an arbitrary function so the general solution of (2.2) is 

q = R b F(cj)-alnR). (2.6) 

Hence in all self-similar disks the pressure, density, velocity of rotation, etc. are 
all of the form (2.6). Notice that S and F are constant on equiangular spirals. 



3. Axially-symmetrical self-similar disks in Equilibrium 
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When we impose axial symmetry the functions F in equation (2.6) all become 
constants. Our equilibrium disks will have a fluid velocity V(R)4> and centrifugal 
force together with the pressure gradient in the disk will balance the gravity of the 
disk. It is useful to define the fictitious circular velocity V C (R)4> as that velocity 
which would by itself balance the gravity — g(R)TL 

+ l f=9(R), (3.1) 

while 

+ ^=9(R) + -^, (3-2) 

R E dR y 1 

where E (-R) is the surface density of the disk and po(R) is the integral of the 
pressure through the (negligible) thickness of the disk. In accordance with equation 
(2.6) we take V c of the form 

V c = B x R- p (3.3) 

and note that because Vq and V c have the same dimensions, Vq is also proportional 
to R~P. Now po/^o has the dimensions of V 2 so it must vary as R~ 2/3 . The 
dimensions of G E are those of [GM/R 2 ] and g, so G E varies as V 2 /R oc .R~( 2/3+1 ). 

Putting these results together 

V 2 /R oc V 2 /R = g oc G E oc R~^+^ , (3.4) 
po = E a 2 oc R-W+V . (3.5) 
Inserting these results into (3.2), we find 

V 2 = V c 2 [l-(4P + l)a 2 /V c 2 ]. (3.6) 
If the pressure decreases outwards a 2 /V 2 < (4(3 + l) -1 . 

Now the gravitational fields of power-law disks are readily calculated e.g. by 
Mestel's method, which gives 

S = 2ttGE /L(/3), (3.7) 

where 

L {0) = T ^~ P)T [\ +f3 } and - \ < (3 < \ . (3.8) 

We have found the approximate formula 

L(J3) = {\- /3)(1.88/3 3 + 1.78/3 2 + 2.44/3 + 2)/(2/3 + 1) (3.9) 

gives a better than 1.2% fit to this function in the required range — \ < f3 < \ 
(Table 1). 
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Table 1 



a 
P 


T ( R\ 
L \P) 


Approximation 


.5 








.4 


.188 


.188 


.3 


.367 


.368 


.2 


.550 


.552 


.1 


.754 


.755 





1.000 


1.000 


-.1 


1.326 


1.329 


-.2 


1.818 


1.829 


-.3 


2.726 


2.755 


-.4 


5.304 


5.348 


-.5 


oo 


oo 



Hence 

2tvG E = L((3)B 2 1 R-^+ 1 ^ = L(p)V?/R (3.10) 

-=^=(S)(^) B " W+1) ' (3 - n) 

For given values of Bi and |/?| < ^ and for a ( constant) chosen ratio of a 2 /V c 2 < 
(4/3+1) > equation (3.3), (3.6), (3.10) and (3.11) define an equilibrium self-similar 
rotating disk. 

The particular interest of the problem 

Notice that the only dimensionful constants defining the equilibrium disks 
are G = [M _1 L 3 T _2 ] and B x = [L /3 + 1 T" 1 ]. The other constants involved (3 
and a/V c are both dimensionless. When we are interested in perturbations about 
the equilibrium we shall take perturbed surface pressures and perturbed surface 
densities to be related by 

Ap/po = 7AS/S0 

where 7 is another dimensionless constant. For j3 7^ 1 it is not possible to make a 
characteristic frequency from the properties of the equilibrium because nothing of 
dimensions [T _1 ] can be made from G and Bi which characterize the equilibrium. 

Now consider an imaginary linear series of equilibria with the same values 
of G, Bi and (3 but with gradually decreasing values of a 2 /V 2 starting from the 
value (4/3 + 1) _1 which is entirely pressure supported. We shall take 7 > 3/2 so 
that the system is then stable. As we decrease the pressure support and increase 
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the rotation speed to compensate, we know we must encounter a system that is 
marginally unstable since the system with a = is unstable to Jeans's instability 
at very short wavelengths. Furthermore, we could perform the experiment of 
exciting modes of only one chosen e tm ^ symmetry. For each m value there will be 
a different marginally unstable mode. Each of these must have a natural frequency 
u m which might or might not be zero. However, u m must be a property of the 
equilibrium configuration and that has no way of making a characteristic frequency 
since no constants with dimensions T _1 can be made from G and B\. We deduce 
that the characteristic frequencies of isolated marginally unstable modes, whatever 
the m may be, must be zero. This implies that the 'pattern' speeds fi p = —u/m 
of isolated modes must also be zero. As we shall see, this greatly simplifies the 
calculations. However, there is another strange possibility originally envisaged by 
Birkhoff in his book. If a mode is not isolated but part of a continuum then it 
can be that the system first becomes unstable with a whole continuum of modes, 
each with a different frequency all becoming unstable together. Such modes would 
be self-similar copies of each other, that of frequency u m being related to that of 
frequency u' m = X~^ +1 ^ui rn by having a scale in space A times as large. In this 
case the system would first go unstable at all scales simultaneously as a result of 
its self similarity. Which of these possibilities occurs? 

For axially-symmetrical modes instability sets in through u = for all invis- 
cid systems with no meridional motions (see e.g. Lynden-Bell & Ostriker). Because 
of this first, Schmitz (1986, 1988 & 1989) and Schmitz & Ebert (1987), and later 
we (Lemos et al. 1991) studied just those modes and found that the stability 
criterion was remarkably close to Toomre's local one. Furthermore, the radial 
displacements of the marginally stable modes were of the form 

in oc R^+V/ 2 cos(a lnR + const) 

or in complex form, £r oc R z where Z = 2(3+ | +i a. For non-axially symmetrical 
modes there is no such nutcracker theorem to indicate that instability should set 
in at u> = but the above arguments indicate that either u = or the mode 
structure is truly remarkable with modes of all real frequencies going unstable 
simultaneously. To spiral structure theorists it would be odd if marginal stability 
took place at u> = for modes of every non-zero m, since this would imply a 
zero pattern speed for the spiral structure; however, the pattern speed is often 
associated with a corotation point at or near the edge of the disk and the power 
law disks are infinite so this could imply a zero pattern speed for them even if this 
were not generally the case. 

These considerations greatly aroused our interest in studying modes with 
frequencies close uj = for the power law disks. As we shall show, the perturbation 
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equations at u> = are self-similar under a radial scaling, so the u> = modes can 
be found exactly for every m. Having found all these modes we then embark on a 
more delicate discussion to discover whether such modes are marginally unstable 
or whether they are ordinary stable modes that happen to have zero frequency in 
non-rotating axes. If the former were true, we would find the marginally unstable 
modes so we could deduce the criterion for stability and the form of the instability 
at each value of m. If the latter were true, we would have proved the remarkable 
result that modes at all frequencies become marginally unstable together so these 
could be sought at any chosen frequency. We shall adopt the Agatha Christie 
approach, not giving the answer until we find it. Those spoil- sports who merely 
want the culprit named can read the end of the paper first. 



4. Disturbances with u> = 

Many years ago now, Lynden-Bell & Ostriker created a powerful formalism 
for discussing stability problems of this type. They showed that the eigenvalue 
problem for u would be reduced to the solution of the equation 

-u 2 A-£ + ooR-£ + Q-£ = (4.1) 

where A, B_ and C_ were Hermitian operators and £e lujt was the displacement 
vector involved in the perturbation. A is positive definite. For the particular case 

in which the undisturbed equilibrium is axially symmetrical with motion in the 4> 
direction only (i.e. no meridional circulation), £ may be taken to vary as e %m ^ and 
the different m components each obey a similar equation with simpler operators 
A , B and C replacing A, B, C. Lynden-Bell and Ostriker's formalism was 
developed for a 3-dimensional system of density po( r ) so we must integrate in the 
vertical direction over all z to obtain equations for a flat disk in which case our 
E (-R) replaces po(r) and all z components vanish. The operator equation then 
takes the simple form 

i-u?A +uB +C )-(l R + %A=0 (4.2) 

where 

4,= E »(i !) < 43 > 

B =-2E„«„( m n +1 ° ,) (4.4) 
=m \ m — 1 J v/ 

where Q (R) = V (R)/R, 

C =T +V +P . (4.5) 

=m =m =m =m v ' 
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Here 



T, = -W< m J 1)J « 1)S | (4.6) 



while and £ m are the operators obtained from y • £ and P • £ by taking £ 
to behave as e im< ^ and making them operate on the vector C = ( ^ ) 

instead of ^ | fl ^ . 

F • £ and P • £ are defined by 

-]/•£ = £ A (WO ( 4 - 7 ) 

£-£ = £ a(^V^ (4.8) 

where ip is the gravitational potential and A is the Lagrangian displacement op- 
erator 

A[Q(r,t)] = Q(T + Z,t)-Q (r,t). (4.9) 

Here Q is any quantity evaluated in the perturbed flow, Q is the same quantity 
evaluated in the unperturbed flow and £ = £(r, t). To first order in £, A is related 
to the Eulerian change 5 evaluated at the same point in both flows by 

A = 5 + £-V. (4.10) 

Lynden-Bell & Ostriker (equation 25) show 

E ■ £ = V [(1 - 7 )Podiv£] - p V(div£) - V [(£ • V)po] + (£ • V) Vp (4.11) 

and from (4.7) and (4.10) 

-Z • £ = S (£ • V)V^ + S V# • (4.12) 

#0 has to be evaluated in terms of £ by use of Poisson's integral remembering that 
5T, = — div(E £). However, although Lynden-Bell & Ostriker used that method 
to demonstrate that V_ was Hermitian, there is a simpler method for the to = 
modes based on their self-similarity. 

The scalars involved in the variational principle of Lynden-Bell & Ostriker 
are quantities such as / / (* ■ T • (RdcfidR. The operator T m itself is proportional 
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to EqO 2 , oc R-^P+V while the RdR integrations raise this to R~^ +1 \ This 
suggests that ( will vary as R+( 4 P+ 1 )/ 2 , so we shall hereafter write 

C = i? (4/3+1)/ V (4-13) 

This agrees with the behaviour found for the axially symmetrical modes by Lemos 
et al. who discuss the interesting behaviour at large R. We write D = Rd/dR 
and after some lengthy calculations done independently by each of us, we obtain 

P ■ C = i? 2/3 - 3 / 2 EoK 2 p -n (4.14) 

— " l =m 

where p • rj = — \o 2 M • ^ and M is the matrix 

/ 7 [ J D 2 -(/3 2 + m) 2 ]+2/3 1 (l + m) 7 [(D - m) 2 - /3 2 ] \ 

^ 7[(D + m) 2_ /3 2 ] 7[jD 2_ (/ 3 2 _ m) 2 ]+2/3i(1 _ m)y ) 

Here 

a 2 = a 2 /y c 2 , /3i = 4/3 + 1 and /3 2 = 2/3 + 3/2 . (4.16) 

Similarly 

T ■ C = R 2l3 ~ 3/2 ^oVnt -rj (4.17) 

where using (3.6) we find 



To calculate the gravity term we see from (4.12) that we need an expression 
for Stf) which arises from 5E. 

SE = -div(E £) = E R- 1 (2/3Cr ~ D^r - im^ R (j>) 

= -E R 2 ^ 1 / 2 [(D + ^r ]R + tm m } ( ' } 



Now E oc iT 2 ^" 1 , so 

R 3/2 5 E oc (D + \) 7] R + irnri^ . (4.20) 

Kalnajs has shown that for flat disks the Fourier transform with respect to u = InR 
of i? 3 / 2 5E is simply related to the Fourier transform in InR of R 1 ' 2 8ip. In particular 

if R 3 / 2 5E oc e lku e im ^, then 

R}/ 2 8*1) = 2kG K(k, m)# 3/2 5E (4.21) 
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where 

Y(m+\+ik)/2)Y{m+\-ik)/2) 
[ ' mJ 5 T((m + f + z/c)/2)r((m + f - ifc)/2) 

For m > 2 Donner has shown this to be well approximated by 

K(k, m) ~ 1/s 

where 



(4.22) 

(4.23) 
(4.24) 



s 2 = (k 2 + m 2 + J). 

To extend Donner's approximation to m = 0, 1 we use the exact recurrence relation 

1 



K(k,m-l)K(k,m) 



(m- \) 2 + k 2 



(4.25) 



which follows from the properties of the Y function. We thus obtain 



K(k, 1) = [K(k, 2) (f + k 2 )] - 1 ~ V /?+¥ (& 2 + |) _1 (4.26) 



tf(M) 1) Q + fc 2 )] 



-l 



fc 2 + ¥ 



fe 2 + | 
^ 2 + l 



(4.27) 



These formulae are accurate to better than 1% as Table 2 shows. Thus when 



rj oc e 



iku-\-imq 



, we find using (4.20) and (4.21) 



Rhi/j = 2ttGKrUj: = -2nGKY, R 2l3 +\ik + \r] R + imr]^} . (4.28) 
(More generally a superposition of components of different ks may be needed) . 

Returning to (4.12) and using (3.10), we finally obtain 

V • f - E V 2 R 2 P =3/2 v • 77 
where (when only one k component is present) 



(4.29) 



v -n = 

=m 



■\LK 



k 2 + (\ + m)' 



(ik — m) 



\-{ik + m) 2 k 2 + (\-m)' 



(3 /3 + 1 
(3+1 (3 



Our operator equation (4.2) can now be written in dimensionless form 



00 u i i 

c c 



v+ 

V- 
(4.30) 



(4.31) 
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where fi c = V c /R, 



" m = [I °) , (4.32) 
Htlo/'m+l , 



and 



a, V m - 1 



c = £ +n + t> . (4.34) 

=m =m — ,™ =m 



Notice that when u = only the c term survives. No term in t , p or v 

=m =m = m =m 

involves u = InR except through the operator D = d/du. Hence these modes have 
an e lku behaviour and we may replace the operator D by ik where it occurs in p . 

=m 

This justifies our use of this form for v for these modes and reduces the problem 
of finding them to algebra. To find such uj = modes we need the determinant 
of to be zero. We find them by calculating the determinant numerically as a 
function of k. Figure 1 plots the k so found against j3. Figure 2 shows the spirals 
that result. Toomre's Q parameter gives the stability of the axially-symmetric 
modes, both in the local approximation (and also quite accurately, < 4%, for the 
global modes of Lemos et al. ). To compare with the Local analysis of a gas disk 
we recall that for such modes 

w 2 = A;2c 2 -27rGEo|A;i| + K 2 (4.35) 

where we use the suffix 1 to distinguish the wave number k\ in R from the wave 
number k in InR used earlier, c is the velocity of sound (7^0/^0)2, k is related to 
fi c via 

k 2 = 40 r 



J__i(Q CJ R 2 ) 
2RdR K ' 



2(1-/3)^. (4.36) 



The expression for u may be rewritten 

.:> _ ,2 



u z =c z \ h 



{ki-^^j + (k 2 — 7t 2 (? 2 Eq/c 2 ) (4.37) 

which demonstrates that the least stable wave number is given by k\ = tt G E / c 2 
and the condition for stability to axially symmetrical modes is Q > 1, where 

kc y/sfl^) c y/8 7 (l - gj 
Q = ^GEo = HP) V c ~ LlJ) a ° • ( ^ 



It is well known that Toomre's original criterion for stars only differs from 
this in that 3.36 replaces the value of 71 in the denominator. Here the gaseous Q 
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defined above is the relevant one. It is also interesting to compare the least stable 
wavenumber k' = koR with the value of k. From (4.35) 



k R = nGZ R/c 2 = \L{f3)V 2 /c 2 = |L(/3)/( 7 ^) 
4(1-/3) 



k R 



(4.39) 



Q 2 Hf3) 



When the modes are non-axially symmetrical, the approximation commonly 
used for tightly wound modes is found from (4.35) or (4.37) by replacing u> by 
uj + mQo(R). Equation (4.35) then takes the form 

~ - ^2m^ + |m 2 (l - (3^1) - [k 2 c 2 - 2vGY> \k\ + ^1 = 

apart from the matrix form of equation (4.31) we see there is a close similarity 
with this approximate dispersion relationship. We therefore expect that uj = 
will occur near where the large curly bracket is zero and that the value of k for 
which uj = will be approximately 



k[ = kiR = 



ro »(i-( V+ i)^-^ + 2Lp> 



Re-writing this expression in terms of our dimensionless variables we find 



k[ 



m\l - (4/3 + 1K)- 2(1-/3) + 



7 2 0- o z l a o 



which is compared with the exactly calculated results in Figure 3. k[ given by 
(4.40) is good approximation for the uj = modes in all but the very open spirals. 
koR shows no agreement. 



5. Marginal Instability 

The existence of stationary modes does not imply marginal stability. To 
ensure the latter we need the limit of growing modes as the growth rate tends to 
zero. Returning to (4.1) we see that unless B_ ■ £ is zero, modes with uj close to 
zero will obey 

uB-i = -c-i 
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However B and C_ are Hermitian and we can see from (4.33) that at least for all 
modes with m > 1, S is negative definite. Hence the eigenvalues uj are necessarily 
real and for any eigen £ 

<£*•£•£> 



u = 



Thus, provided < £* • 5 • £ >^ 0, all the eigen frequencies close to u = are 
real. Now we may evaluate < £* ■ B_ ■ £ > for our uj = mode it can only 
be slightly different for the infinitesimally different mode with uj — > 0. Hence, 
provided < £* • S • £ >^ for our u; = modes, there can be no marginally 
unstable modes close to ui = . 

Now in the form (4.33) b m is diagonal, so all we need to ensure is that for 
our u = mode 

multiplying by C22 the LHS becomes 

V+V+ (&11C22 + &22C11) 

where we have used the fact that c m ■ rj = to change the r?_ terms to r? + ones by 
writing 

C22V- = - c i2V+ anci c 12 r)- = -cnr/+ . 

Thus, proveded 77+ is non zero, the condition 611C22 + &22C11 7^ implies that 
< £, ' B.' £, >7^ f° r the marginal mode. (Even if 77+ is zero, a similar argument 
multiplying by en to start with yields the same result provided rj- is non zero and 
they can not both be zero as then 77 vanishes). Now 



O11C22 + O22C11 = — - — 



-2(m 2 - 1 + 0) - L K(m 2 + k 2 - |) + 
+a 2 {2/3 1 (m 2 - 1) + 7 [(/c - /3 2 ) 2 + m 2 ]} 



To find out whether this can ever be zero for a marginally stable mode we solved 
the equation 611C22 + &22C11 = for a 2 and inserted the result into ||c m || for each 
k, the result was negative for every k and every m^O. 

We deduce that any marginally stable mode with m 7^ has a non-zero 
frequency ui m and an associated corotation resonance (provided ui m is positive) 
at = B\/u} m . The re-scalings of this mode by factors A yield an infinity 

of marginally stable modes with frequencies X~^ +1 'uj m . Thus, marginally stable 
modes occur simultaneously at all frequencies. We find this result to be unexpected 
and remarkable. One may now search for a marginally unstable mode by first 
choosing a frequency; any one is as good as any other; it is unnecessary to search 
through all possible frequency as one would have to do for a non-self-similar system. 
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Axially Symmetric Modes 

It is of interest to see how the uj = modes studied earlier fit into the 
present analysis. Setting m = 0, c becomes a multiple of the matrix 



1 1 

and 611C22 + 622C11 vanishes. Thus the determinant \\c m \ \ and < • S • £ > both 
vanish. This is associated with the existence of the trivial neutral displacement 
with £^ independent of 0. To find the marginally unstable modes one must expand 
the determinant 

f 00 \ u , 

- — a + 7^6 +c 



to all orders in uj/Q c when one recovers 

w 4 ,, M f u \ 3 ( uj \ 2 f anc 2 2 + a 2 2Cn 

7^7 m - FT ( a ii»22 + a220n) - , , 

but for m = 0, 011622 + ^22611 = 0, so the equation for the non-trivial uj = 
axially symmetric modes is 

ui 2 

^2l|a m || = (anc 2 2 + a22Cn - 611622) = 

For m = 0, anc 2 2 + fl22Cn - 611622 = yields 

°l\l(k 2 + (3D - 4/3i] + 2(1 -0)- LK(k 2 + \) = 

which is the condition of marginal stability for m = modes found by Lemos et 
al. in a slightly different notation. 

Conclusion. 



In line with the anti-spiral theorem of Lynden-Bell and Ostriker the condition 
||c m || = involves k 2 , m 2 and VL 2 but not the signs of k, m, fi c . This implies that 
to every trailing spiral mode with u> = there corresponds a leading spiral mode 
with uj = and that we could take pairs of linear combinations of them which 
are not spiral at all! Only for growing or decaying modes can this symmetry be 
broken in a dissipationless system. 

The non-axially symmetric uj = modes are not the key to the stability 
problem, rather they are examples of stable steady modes which can be found 
exactly. However, that result leads to the remarkable conclusion that modes of all 
frequencies become marginally unstable together, so that marginally stable modes 
can be sought at any frequency and others can then be found by self-similar scaling. 
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Table 2 



1 

K 


K(k,0) 


Approx 


TV" / 1 -1 \ 

K(k,l) 


Approx 


K(k,2) 


Approx 





4.377 


4.366 


0.914 


0.916 


0.486 


0.485 


0.2 


3.834 


3.812 


0.9018 


0.904 


0.484 


0.483 


0.4 


2.810 


2.799 


0.868 


0.871 


0.478 


0.476 


0.6 


2.004 


1.993 


0.818 


0.823 


0.468 


0.466 


0.8 


1.479 


1.468 
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Figure Captions 



Figure 1. k the dimensionless wavenumber in InR space for u> = waves in 
power law disks with rotations V oc . The 'ring' modes m = are marginally 
unstable. The m^0,w = modes are stable. 

Figure 2a. u = modes of V oc R -13 gas disks: On the left (3 = 0, Q = 1, 7 = 2; 
bottom m = 1, /c = 7; second m = 2, = 9.45; third m = 3, /c = 11.9; top m = 4, 
= 14.5. On the right /? = 0, Q = 1.5, 7 = 2; bottom m = 1, k = 2.6; second 
m = 2, /c = 4.6; third m = 3, k = 6.1; top m = 4, k = 7.7. Notice that the larger 
<5 gives more open spirals. 

Figure 2b. (3 = —0.2. For equal Q and m the less sheared rotation law gives 
more open spirals Q = 1, 7 = 2, m = 2, A; = 5.8. 

Figure 3. The dimensionless wavenumber of the exact u; = mode is plotted 
against k' , the value of koR for marginal stability in the local tightly- wound ap- 
proximation, k and k' give disparate curves k[, the fci-R of the uj = stable mode 
in the tightly wound approximation is a reasonably good approximation for large 
k. 
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This figure "2a.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9907093vl 



This figure "2b.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9907093vl 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9907093vl 



This figure "fig3.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9907093vl 



